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The four wave interaction process, known as the oscil- 
lating two stream instability (OTSI) is considered as one 
of the mechanisms responsible for stabilizing the electron 
beams associated with solar type III radio bursts. It has 
been reported that (1) an intense localized Langmuir wave 
packet associated with a type III burst contains the spec- 
tral characteristics of the OTSI: (a) a resonant peak at the 
local electron plasma frequency, f pe , (b) a Stokes peak at a 
frequency slightly lower than / pe , (c) anti-Stokes peak at a 
frequency slightly higher than / pe , and (d) a low frequency 
enhancement below a few hundred Hz, (2) the frequencies 
and wave numbers of these spectral components satisfy the 
resonance conditions of the OTSI, and (3) the peak intensity 
of the wave packet is well above the thresholds for the OTSI 
as well as spatial collapse of envelope solitons. Here, for 
the first time, applying the trispectral analysis on this wave 
packet, we show that the tricoherence, which measures the 
degree of coherent four-wave coupling amongst the observed 
spectral components exhibits a peak. This provides an ad- 
ditional evidence for the OTSI and related spatial collapse 
of Langmuir envelope solitons in type III burst sources. 


1. Introduction 

Solar type III radio bursts are characterized by very fast 
negative frequency drifts from hundreds of MHz to tens of 
kHz. Ginzburg and Zheleznyakov [1958] were the first to 
propose that the production of these bursts involves the 
excitation of high levels of Langmuir waves at local elec- 
tron plasma frequency f pe = 971*, where n e is the electron 
density in m -3 by flare accelerated electron beam through 
bump-on-tail instability [Bohm and Gross , 1949], and their 
subsequent conversion into radio emissions at fpe and 2/pc , 
which has been confirmed by the in situ detection of Lang- 
muir waves [Gumett and Anderson , 1976, 1977] as well as 
electron beams [Lin, 1970; Lin et al. , 1973, 1986]. According 
to Sturrock [1964], the excitation of Langmuir waves would 
extract all the streaming energy from the electron beam 
within 100 km or less, whereas, the bump-on-tail distribu- 
tions of electron beams are detected [Lin et al . , 1986] over 
the distances of 1 AU and more. It is now believed that some 
nonlinear process removes the Langmuir waves rapidly from 
the spectral regions of resonance with the beam, which leads 
to the beam stabilization. For example, the induced scat- 
tering off ion clouds , which is the electrostatic decay (ESD) 
of initial Langmuir wave into a daughter Langmuir wave 
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and an ion sound wave when T e > T* [Bardwell and Gold- 
man, 1976] is proposed as one of such mechanisms [Kaplan 
and Tsytovich , 1968], where, T e and Ti are the electron and 
ion temperatures, respectively. Although the signatures of 
electrostatic decay are observed in the type III sources [Lin 
et al ., 1986; Gumett et al ., 1993; Hospodarsky and Gumett , 
1995; Thejappa and MacDowall , 1998; Thejappa et al. , 2003; 
Henri et al., 2009], their time scale appears to be too long 
to prevent the plateau formation [Zheleznyakov and Zaitsev , 
1970). 

The type III associated Langmuir waves are usually esti- 
mated to be very intense and therefore, the four-wave inter- 
action called the oscillating two stream instability (OTSI) 
[Papadopoulos et al., 1974; Smith et al., 1979; Goldstein et 
al., 1979], and related soliton formation and spatial col- 
lapse [ Zakharov , 1972; Nicholson et al., 1978] are proposed 
as the most effective beam stabilization mechanisms. The 
OTSI excites a low frequency ion density perturbation of 
frequency and wave number ($7, q), which can beat with 
two of the beam-excited Langmuir waves of frequency and 
wave number ( f pe , /cl) and produce down-shifted ( f pe — f2, 
kj, — q) (Stokes) and up-shifted (/ pe + Q, hL+q ) (anti-Stokes) 
modes, respectively. The spatial collapse, on the other hand 
occurs due to intensification of the localized Langmuir wave 
packet in the self generated shrinking density cavity. Some 
possible evidence for the strong turbulence processes in the 
Jupiter’s foreshock [Gumett et al., 1981], in the solar wind 
[Kellogg et al., 1992], and in the source regions of type III 
bursts [Thejappa et al., 1993; Thejappa and MacDowall , 
1998; Thejappa et al., 1999; Thejappa and MacDowall , 2004] 
was reported. 

In a recent study, Thejappa et al. [2012] have reported 
the STEREO/S WAVES [Bougeret et al., 2008] high time 
resolution observations of an isolated localized type III as- 
sociated Langmuir wave packet with short duration of ~ 3.2 
ms. These authors have shown that (1) the spectrum of this 
wave packet contains the characteristics of OTSI, namely, an 
intense peak at f pe , and two side bands at slightly lower and 
higher than f pe , and a low frequency enhancement below a 
few hundred Hz, (2) the frequencies and wave numbers of 
these spectral components satisfy the resonance conditions 
of the OTSI, and (3) the peak intensity of the wave packet 
is well above the thresholds for the OTSI as well as for the 
formation of envelope solitons collapsed to a few hundred 
Debye lengths. Based on these observations, it has been ar- 
gued that the OTSI and spatial collapse control the beam 
plasma interactions in the solar type III radio bursts. 

In this study, for the first time, we will apply trispec- 
tral analysis on this Langmuir wave packet and compute 
the tricoherence. We will show that the tricoherence spec- 
trum exhibits a peak, which is an indicative of the four-wave 
interaction of type OTSI. We argue that the high degree 
of phase coherence between the spectral components of the 
wave packet provides an additional evidence for the OTSI 
and related strong turbulence processes in the solar type 
III radio bursts as correctly concluded by Thejappa et al. 
[2012]. In section 2, we review the observations, in section 
3, we present the trispectral analysis and in section 4, we 
present the conclusions. 
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2. Review of Observations 

In Fig. 1, the fast drifting emission from very high fre- 
quencies to the local electron plasma frequency, f pe ~ 30 
kHz is identified as the local type III burst, and the non- 
drifting emissions in the interval 27-32 kHz are identified as 
the Langmuir waves. Time Domain Sampler (TDS) of the 
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Figure 1 . Dynamic spectrum of a local type III radio 
burst (fast drifting emission from ~ 5 MHz down to ~ 30 
kHz) and associated Langmuir waves (non-drifting emis- 
sions in the frequency interval 27-32 kHz). 


STEREO A - 2010:09:12:07:08:18.342 



0 10 20 30 40 50 60 

Time (ms) 


Power Spectrum 



Frequency (kHz) 


Power Spectrum 



Frequency (kHz) 


Figure 2. (a) The Langmuir wave packet observed 

by the Time Domain Sampler (TDS) during the type 
III event of Fig. 1, (b) The narrow spectrum around 
/ ~ fpe 30 kHz. The L, S, and AS correspond to the 
beam excited Langmuir wave, Stokes peak at ~ 29.54 
kHz, and anti-Stokes peak at ~ 30.41 kHz, respectively, 
and (c) Low frequency spectrum: the enhancement below 
450 Hz is probaby due to ion-sound waves 


S WAVES experiment [Kellogg et al. , 2009], which samples 
the A/C electric field from 3 orthogonal antennas has re- 
solved these Langmuir waves into intense waveforms, each 
of which contains 16384 samples with an acquisition rate of 
250,000 samples per second (a time step of 4 /jls for a total 
duration of 65 ms). The most intense wave packet, captured 
by the E x antenna is shown in Fig. 2a. The peak electric 
field strength El and ^-power duration r of this event are 
56.5 mVm" 1 and ~ 3.2 ms, respectively. The narrow spec- 
trum of this event, as seen in Fig. 2b shows an intense peak 
(L) at fpe ~ 30 kHz, corresponding to n e ~ 1.1 x 10 7 m“ 3 , 
a Stokes peak (S) at ~ 29.54 kHz, which is slightly less than 
fpe and an anti-Stokes peak (AS) at ~ 30.41 kHz, which 
is slightly higher than f pe . The spectrum, as shown in Fig. 
2c clearly shows a low frequency enhancement corresponding 
to ion-sound waves below 450 Hz. The STEREO/PLASTIC 
experiment [Galvin et al, 2008] has measured the solar wind 
speed v S w as ~ 450 kms' 1 . We assume that the electron 
temperature T e is ~ 10 5 K during this event. Assuming that 
the type III electrons propagate along the Parker’s spiral 
field lines, we fit a frequency drift curve to the dynamic 
spectrum and estimate the beam speed Vb as ~ 0.22c for 
the RAE density model [Fainberg and Stone , 1971], where 
c is the velocity of light. However, the pitch angle scatter- 
ing is known to increase the path length of electron beams 
by a factor of a = 1.3 to 1.7 [Alvarez et al, 1975; Lin ct 
al , 1973]. This implies that if we incorporate these cor- 
rections into the estimates of the beam speeds, they will 
lie in the range from ~ 0.29c to ~ 0.37c. Accordingly, the 
wave number of the Langmuir waves /cl = ^ ~ 2.9 x 10 -3 
m -1 estimated for Vb = 0.22c will decrease to the values 
from ~ 2.2 x 10 -3 m“ x to ~ 1.7 x 10 -3 m _1 . We es- 
timate the rest of the parameters as: (1) Debye length, 
A De = 69Te ,2 n7 l/2 ~ 6.6 m, and (2) the normalized peak 
energy density — 10“ 3 , and (3) the —power 

spatial scale of the wave packet S ~ 219Aoe (using the re- 
lation S ~ TVsw)- 

It was argued that these observations of a strong Lang- 
muir wave peak with upper and lower sidebands (Fig. 2b), 
together with low frequency enhancement (Fig. 2c) are 
strongly suggestive of OTSI, in which, the beam driven 
Langmuir v. ave is the pump wave, the modes corresponding 
to sidebands and low frequency waves are the nonlinearly 
excited daughter waves. The frequency matching condition 
2 /l =/s + /as is easily satisfied, since the frequency shifts 
of the Stokes and anti-Stokes modes are symmetric with re- 
spect to the Langmuir wave pump, being ~ 442.5 Hz and 
~ 427 Hz, respectively, which in turn are in good agreement 
with the frequencies of the observed ion sound waves of < 
450 Hz. The wave numbers k of the side bands are estimated 
using the expression for the frequency shift [Gurnett et al , 
1981] A/= 2^(/cA De )cos6»-F/ pe (-l-h(l+3(/cADe) 2 ) 1/2 ), 
where, the first and second terms correspond to the Doppler 
shift caused by the motion of the solar wind and intrinsic 
frequency variation caused by the dispersion relation, re- 
spectively, and 0 is the angle between k and 0 = 0 
and 0 = 7T correspond to the anti-Stokes and Stokes modes 
propagating away from and toward the Sun, respectively. 
By plugging the measured A / of ~ 442.5 Hz and 427 Hz 
in this equation, kXp is estimated as ~ 0.03 and ~ —0.05, 
for the anti-Stokes and Stokes modes, respectively, which 
suggest that the pump Langmuir waves with fc^Aoe ~ 10 -2 
are pumped into those of forward and backward propagat- 
ing daughter waves with large wave numbers. The upper 
limit of the wave numbers of the ion sound waves can be 
estimated using the relation q = ^ ZL ^, since their phase ve- 
locities are usually well below u stjU f u Thus, for Q = 450 Hz 
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and v 3W = 450 kms -1 , it is estimated that q ~ 5.6 x 10 -3 
m _1 and qXoe — 0.036. These wavenumbers are compara- 
ble to those of the sideband emissions, and the matching 
condition fc = /cl ± g is easily satisfied, yielding \k\ ^ |g|, 
since /cl << <?• 

The threshold for the OTSI [ Zakharov , 1972] > 

{kL^De) 2 is easily satisfied, since the observed ~ 10” 3 
is well above (/clAdc ) 2 ~ xl0“ 4 . For the wave packet 
to be the collapsed soliton, it should satisfy the condi- 
tion [Thornhill and ter Haar, 1978; Gwnett et al. , 1981] 
> (AkXpe) 2 , where A k = is the wavenumber char- 
acteristic of the envelope. In the present case, the observed 
~£- ~ 10 -3 is greater than (A/cAdc) 2 ~ 8 x 10” 4 obtained 
for the spatial scale S ~ 219Ad c . This suggests that the 
observed wave packet is probably the Langmuir envelope 
soliton, collapsed to the spatial scale of ~ 219 A£> e - 

3. Trispectral Analysis 

The phase coherence between the spectral components is 
one of the important characteristics of the four- wave interac- 
tion, such as the OTSI. In the power spectrum estimation, 
the waveform is treated as a superposition of statistically 
uncorrelated harmonic components, and the phase relations 
between the spectral components are suppressed. The in- 
formation present in the power spectrum is sufficient for the 
complete statistical description of any Gaussian process of 
a known mean. However, in order to extract the informa- 
tion regarding the presence of nonlinearities, we should look 
beyond the power spectrum. Higher order spectra (HOS), 
which are defined in terms of the higher order moments or 
cumulants of the signal contain such information. The third- 
order spectrum is commonly referred to as bispectrum [ Kim 
and Powers , 1979], the fourth-order one as trispectrum. The 
information about the phase coherence between four spectral 
components can be extracted from the waveform data using 
the trispectral analysis. The trispectral method has been 
developed and applied to synthetic [Kravtchenko -Berejnoi 
et al., 1995a, b; Lefeure et al., 1995] and to simulated data 
[Soucek et al., 2003]. The trispectral analysis can detect the 
phase relationships among four Fourier components, which 
is the key information regarding four-wave interactions. In 
OTSI, the sidebands interact with the strong beam excited 
Langmuir waves simultaneously satifying the matching rules 
for the wavenumbers as well as for frequencies. For the OTSI 
type of four-wave interactions, the cumulant based trispec- 
trum is given by [Kravtchenko -Berejnoi et al., 1995a] 


T(l, 2, 3) = E[X l X 2 X;x:\ - N( 1,2, 3, 4), (1) 

where (X\,X 2 ,Xz) and X 4 are the complex Fourier 
components of the signal corresponding to frequencies 
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Figure 3. The cross-section at / = 29.5 kHz of the tri- 
coherence spectrum of the TDS event of Fig. 2a. 


fu h, h and / 4 , TV (1, 2, 3,41 = E[X l X 2 ]E[X^X:] + 
E^X^E^XZ] + E[XxXZ)e[x 2 XZ), ft = f x + f 2 - / 3 , 
and E\\ is the expectation operator. The tricoherence, which 
is the normalized trispectrum is usually used, because it 
eliminates the dependence of trispectrum on the ampli- 
tude of the signals. The expressions for cumulant based 
square tricoherence function can be written as [Kravtchenko- 
Berejnoi et al., 1995a]: 

t 2 n 2 3 ) = l r d.2,3)| 2 (2) 

M,,; (£[|XiX 2 X 3 *X 4 *|] + |W(1,2,3,4)|) 2 ‘ {l) 

A unit value for the tricoherence indicates perfect coupling, 
a zero value indicates no coupling, and any \alue between 
zero and one indicates partial coupling. The triphase is the 
phase of the tricoherence. The tricoherence quantifies the 
fraction of the total product of powers at the frequency quar- 
tet, (/1, f 2 , fs 1 f\ + /2 - / 3 ), that is owing to cubicly phase- 
coupled modes. The tricoherence is zero for a Gaussian pro- 
cess. 

For Gaussian process, due to statistical fluctuations, the 
estimate of the tricoherence from a finite data record will 
not be zero. Therefore, the method of periodograms is usu- 
ally used to estimate the trispectrum and tricoherence. This 
involves the division of the data record of N samples into M 
time intervals and calculation of the ensemble averaged es- 
timates over all intervals. In this study, we have used the 
segment length N = 1000 (0.004s), number of segments M 
= 16, Hamming windowing, and a cumulant estimator. We 
have calculated the tricoherence spectrum as a function of 
three frequencies. Since it is difficult to visualize the results 
in such a 3-D space, we have made the cross-section of the 
tricoherence domain at the frequency / 3 = 29.5 kHz, corre- 
sponding to the frequency of the Stokes mode. In Fig. 3, we 
present such a cross section. The tricoherence peak shown 
as circled feature at (/1 = 30 kHz, f 2 = 30 kHz, / 3 = 29.5 
kHz) is due to the phase relation 2^1, = <j>s + where 
4>s and (f>AS are the phases of the beam-excited Lang- 
muir, Stokes mode and anti-Stokes modes, respectively. The 
maximum tricoherence t 2 is ~ 0.51. This significant value 
of tricoherence provides evidence for the cubic nonlinearity, 
specifically for four- wave interaction, known as the OTSI. 

4. Conclusions 

Thejappa et al. [2012] reported that (1) TDS of the 
SWAVES experiment on STEREO A has captured very co- 
herent and intense Langmuir wave packets in the source re- 
gion of a local type III radio burst, (2) The spectrum of this 
wave packet contains the characteristic signatures of oscil- 
lating two stream instability (OTSI), namely, (a) a resonant 
peak at the local electron plasma frequency, / pe , (b) Stokes 
peak at a frequency slightly lower than f pe , (c) anti-Stokes 
peak at a frequency slightly higher than /pe> and (d) low 
frequency enhancement corresponding to ion sound fluctu- 
ations; the frequencies and wave numbers of these spectral 
components satisfy the resonance conditions of OTSI, (2) 
The observed ~ 10“ 3 , which is well above the thresh- 
old for strong turbulence processes, and the short time scale 
~ 3.2 ms indicate that the observed wave packet is a collaps- 
ing Langmuir envelope soliton, (3) the Langmuir collapse 
follow’s the plane-wave modulational instability as studied 
by Zakharov [1972], and (4) the strong turbulence processes 
are probably responsible for beam stabilization as well as 
for the observed type III burst fundamental and harmonic 
emissions. 

In this study, for the first time, using the trispectral anal- 
ysis, we have shown that the tricoherence between the beam 
excited Langmuir wave and the side bands is ~ 0.51, which 
provides evidence for the four wave interaction of the type 
OTSI and provides an additional support for the conclusions 
drawn by Thejappa et al. [2012]. 
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